Compensated Isocurvature Perturbations and the Cosmic Microwave Background 
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Measurements of cosmic microwave background (CMB) anisotropics constrain isocurvature fluc- 
tuations between photons and nonrelativistic particles to be subdominant to adiabatic fluctuations. 
Perturbations in the relative number densities of baryons and dark matter, however, are surprisingly 
poorly constrained. In fact, baryon-density perturbations of fairly large amplitude may exist if they 
are compensated by dark-matter perturbations, so that the total density remains unchanged. These 
compensated isocurvature perturbations (CIPs) leave no imprint on the CMB at observable scales, 
at linear order. B modes in the CMB polarization are generated at reionization through the mod- 
ulation of the optical depth by CIPs, but this induced polarization is small. The strongest known 
constraint < 10% to the CIP amplitude comes from galaxy-cluster baryon fractions. Here it is shown 
that modulation of the baryon density by CIPs at and before the decoupling of Thomson scattering 
at 2 ~ 1100 gives rise to CMB effects several orders of magnitude larger than those considered 
before. Polarization B modes are induced, as are correlations between temperature/polarization 
spherical-harmonic coefficients of different Im. It is shown that the CIP field at the surface of last 
scatter can be measured with these off-diagonal correlations. The sensitivity of ongoing and future 
experiments to these fluctuations is estimated. Data from the WMAP, ACT, SPT, and Spider 
experiments will be sensitive to fluctuations with amplitude ~ 5 — 10%. The Planck satellite and 
Polarbear experiment will be sensitive to fluctuations with amplitude ~ 3%. SPTPol, ACTPol, 
and future space-based polarization methods will probe amplitudes as low as 0.4% — 0.6%. In the 
cosmic- variance limit, the smallest CIPs that could be detected with the CMB are of amplitude 
~ 0.05%. 

PACS numbers: 98.70.Vc,95.35.-|-d,98.80.Cq,98.80.-k 



I. INTRODUCTION 



The concordance cosmological model posits a nearly 
scale-invariant spectrum of primordial density fluctua- 
tions with adiabatic initial conditions, for which the ra- 
tios of neutrino, photon, baryon, and dark-matter num- 
ber densities are homogeneous. The simplest inflationary 
models predict adiabatic fluctuations ^^ , and adiabatic 
fluctuations are consistent with measurements of cosmic 
microwave background (CMB) temperature/polarization 
anisotropies [3] and the clustering of galaxies [1, Q ■ 

Isocurvature perturbations are fluctuations in the ra- 
tios of number densities of various particle species. They 
are produced in topological-defect models for structure 
formation UM and in more complicated models of infla- 
tion [ll| - [lq |. CMB temperature anisotropies limit the 
amplitude of baryon isocurvature perturbations (fluctua- 
tions in the baryon-to-photon ratio) [13, \l3, and CDM 
isocurvature perturbations (fluctuations in the dark- 
matter-to-photon ratio) Jl*[22[to be < 13% of the total 
perturbation amplitude [7|, I23l433l | . 



turbations (CIPs) thus obey 



CI 



pj^' + pb<5^^ - 



xci 



-0, 



(1) 



Our intuition thus suggests the matter in the early Uni- 
verse was very smoothly distributed. It therefore comes 
as somewhat of a surprise to learn that perturbations 
in the baryon density can be almost arbitrarily large, as 
long they are compensated by dark-matter perturbations 
such that the total nonrelativistic matter density remains 
unchanged [3J,|35|. These compensated isocurvature per- 



where Sc, Sb, and S^ are fractional energy density per- 
turbations in the dark matter, baryons, and photons, re- 
spectively, while pc and pb are the homogeneous dark 
matter and baryon densities. 

CIPs induce no curvature perturbation at early times, 
and they therefore leave the photon density — and thus 
large-angle CMB fluctuations — unchanged at linear or- 
der. CIPs induce baryon motion through baryon-pressure 
gradients, but these motions occur only at the baryon 
sound speed which, at the time when Thomson scatter- 
ing first decouples (z ^ 1100, decoupling hereafter), is 
{v/c) ~ (T/mp)i/2 - (eV/GeV)i/2 ^ 10-4-5. The ef- 
fects of these motions on CMB temperature and polar- 
ization anisotropies thus occur only on distances smaller 
than ~ IQ-^-S times the sound horizon at decoupling or 
CMB multipole moments / ~ 10^ [13, [H, 113, scales far 
smaller than those probed by CMB experiments. 

The effect of CIPs on galaxy surveys is also believed 
to be small [33|. Big-bang nucleosynthesis (BBN) and 
galaxy-cluster baryon fractions constrain the CIP per- 
turbation amplitude to be < 10% [3a|. Measurements 
of fluctuations in 21-cm radiation from atomic hydrogen 
during the dark ages may be sensitive to these perturba- 
tions [34, 37 39], but these measurements are a long way 
in the future. 

In Ref. [33, it was shown that, although CIPs produce 



no observable effect on the CMB at linear order in per- 
turbation theory, they modulate the CMB fluctuations 
produced by adiabatie perturbations. In particular, it 
was shown that B modes in the CMB polarization are 
produced by the angular modulation in the reionization 
optical depth induced by the CIP. 

Here, we consider the additional effects on CMB fluc- 
tuations that arise from modulation of the baryon den- 
sity by CIPs at and before decoupling. CIPs modulate 
the free-electron density. They thus change the photon 
diffusion length and thickness of the surface of last scat- 
tering (SLS) on different patches of sky. CIPs also change 
the weight of the baryon-photon plasma and thus the de- 
tails of the acoustic-peak structure in the CMB power 
spectrum. Variation in the baryon density from one re- 
gion on the sky to another thus leads to a modulation 
of the small-scale power spectrum from one region of 
sky to another. This induces B modes in the polariza- 
tion and nontrivial higher-order correlations in the tem- 
perature/polarization map analogous to those induced 
by variations of other cosmological parameters (40[ and 
those induced by weak gravitational lensing [4l| . 

As we show below, the effects of CIPs on CMB fluctu- 
ations from decoupling are several orders of magnitude 
larger than those from reionization, and so the CMB 
should provide a far more sensitive probe of CIPs than 
envisioned in Ref. [35|- We therefore follow through 
and develop the formalism required to look for CIPs 
with the CMB. To do so, we write down the minimum- 
variance estimators that can be constructed from a CMB 
temperature-polarization map for the CIP field A(n) as 
a function of position h on the sky. We evaluate the 
noise with which the CIP field can be reconstructed and 
estimate the signal-to-noise with which a scale-invariant 
spectrum of CIPs may be detected with various experi- 
ments. 

We conclude that data from WMAP, Spider, ACT, and 
SPT are sensitive to CIP amplitudes of - 5 - 10%. The 
Planck satellite |42| and Polarbear experiment are sen- 
sitive to CIP amplitudes as small as ~ 3%. Upcoming 
ground-based polarization experiments (ACTPol [ij] and 
SPTPol [Mlill) or a post-Planck CMB-polarization ex- 
periment along the lines of the proposed EPIC experi- 
ment [Hi could detect fluctuations of ^ 0.4% - 0.6%. In 
the cosmic-variance limit, sensitivity to fluctuations of 
amplitude ^ 0.05% is possible. 

Our principal motivation in studying CIPs is curios- 
ity: can we determine empirically, rather simply assume, 
that the primordial baryon fraction is homogeneous and 
traces the dark matter? Still, there may be theoreti- 
cal motivation as well. For example, curvaton models 
for inflation may generate CIPs [Sa, |47H49| , with ampli- 
tudes approaching the regime detectable by EPIC [M] • It 
may also be that recent models |50l - [55| that connect the 
baryon asymmetry and dark-matter density have impli- 
cations for CIPs. Additionally, the techniques introduced 
in this paper could be used to empirically disentangle 
a CDM isocurvature fluctuation from a baryon isocur- 



vature fluctuation, using CMB data. These modes are 
usually treated as degenerate in the analysis of CMB ob- 
servations. 

In Ref. [56|, we presented our basic conclusions. Here 
we present in detail our results, their derivation, and the 
computational methods used. We calculate the induced 
temperature anisotropics in Sec. |TT] and the induced po- 
larization anisotropies in Sec. IIIIl In Sec. IIVI we com- 
pute the expected corrections to CMB power spectra for 
a scale-invariant spectrum of CIPs and compare the B- 
mode power spectrum induced by CIPs at decoupling 
with that induced at reionization. In Sec. |Vl we con- 
struct minimum-variance estimators for CIPs. We then 
assess in Sec. I VII the sensitivity of ongoing and upcom- 
ing experiments to CIPs, and we conclude in Sec. I VIII 
Useful relations involving tensor spherical harmonics are 
presented in Appendix 1X1 Numerical derivatives of trans- 
fer functions arc discussed in Appendix [BJ Second-order 
harmonic expansions for CMB transfer functions are de- 
rived in Appendix [C] Throughout, we use as our fiducial 
cosmological parameters those from Ref. [3] . 



II. PERTURBED LINE-OF-SIGHT (LOS) 
FORMALISM: TEMPERATURE 



Here we review the standard calculation of the 
temperature-fluctuation power spectrum for primordial 
adiabatie density perturbations. We then show how this 
calculation is altered in the presence of CIPs. 



A. General line-of-sight solution for temperature 

The spherical-harmonic coefficients Tim for the CMB 
temperature can be written 



Tim = / dnT{n)Yim{n) 

dnyi*m{n)Yi^m,{.n) / dr]f{r],n) 
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where T (n) is the CMB temperature in direction n, and 
Yimin) are spherical harmonics. The Fourier transform 
of the primordial gravitational potential for wave-vector 
k is $^, while ji{x) denotes a spherical Bessel function. 
The conformal time rj = J dt/a(t) is here an integration 
variable, and ryo denotes its value today. The function 
/ (77, n), obtained via the numerical solution of the Boltz- 
mann equations |57H59l |. encodes how much a real-space 
primordial-potential perturbation $[(77 — 770)71,77] con- 
tributes to the temperature anisotropy T{n). It depends 
on the relation between initial gravitational-potential 
fluctuations and radiation-density fluctuations at decou- 
pling, as well as the recombination history. 



B. Temperature anisotropies with homogeneous 
baryon fraction 

In the standard calculation, this transfer function is 
the same in all directions; i.e., f(r],h) = /(t?). In this 
case. Eg. ([2l) simplifies, via orthogonality of the YimS, 
yielding [58|, |60| 



Tim — 



47ri' 
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dvfiv) / d'k'^j:ji[k{m-v)]YCUk). (3) 



The temperature power spectrum is then easily obtained 
using Eq. (J3|), averaging over realizations of the poten- 
tial perturbation, and using the identity ($(fc)$*(fc')) = 

{27Tf S^jj{k-k')P.s>{k), where 5l,{k-k') is the Dirac delta 
function and P^{k) the primordial-potential power spec- 
trum, and the angle brackets denote an average over re- 
alizations of the primordial potential. We then find [S^] 



{Tl'm'T^lm) - Ci Su'S, 



where 
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(5) 



is the CMB temperature power spectrum, written in 
terms of a transfer function. 



Ti{k) 



dvfiv)Ji[Hvo-v)]- 
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This transfer function is tabulated by Boltzmann codes 
like CAMB [m and CMBFAST [sll , and Sij is the Kronecker 
delta. 



Temperature anisotropies with CIPs: Single 
CIP realization 



In the presence of a compensated isocurvature pertur- 
bation, the baryon and dark-matter fractions vary from 
one point in the Universe to another, and so the transfer 
function f{rj,n) now acquires some direction (n) depen- 
dence. The CIP involves small changes. 



rib ^ rib [1 + A (n)] , 
r^c — >■ ric - fibA (n) , 



(7) 



in the cosmological parameters between different lines of 
sight h. Here, A{n) is the value of the CIP in direction 
n at the surface of last scatter (or reionization — wc will 
make these statements more precise below). Note that 
we define it so that it is the fractional perturbation in 
the baryon (rather than dark-matter) density associated 
with the CIP. From Eq. ([7]) , the change in the total den- 
sity is (5 rim = i^^c + S^h = 0, and so this is indeed a 
compensated isocurvature perturbation. 

In a general treatment of perturbed recombina- 
tion/decoupling, one would follow the set of equations 



for the electron, dark-matter, photon, and neutrino den- 
sities, velocities, and t he g ravitational potential at second 
order, as in Refs. [63 - [64| . In the case of CIPs, however, 
the CIP amplitude does not evolve for all observation- 
ally accessible scales, and we can thus model the effect 
of CIPs as a modulation in the cosmological parameters 
ric and r^b. 

We can build some intuition for the effect of CIP per- 
turbations on the CIP by considering a globally constant 
CMB perturbation A. We run the CAMB [6l| code with 
a global perturbation of the form in Eq. (O for a variety 
of A values. We evaluate the angular sound horizon Ig at 
the surface of last scatter as a function of A, using the ex- 
pressions in Ref. [631 ■ W^e see in the top left panel of Fig. 
[T]that, as the plasma is more loaded down with baryons 
in the presence of a CIP with a positive A value, the de- 
crease in sound speed moves the CMB acoustic peaks to 
smaller angular scales. 

CMB temperature anisotropies are suppressed on an- 
gular scales / > Id due to diffusion damping. Using the 
expressions in Ref. [60| and the CAMB 61 1 code, we eval- 
uate Zd(A) and show the results in the top right panel 
of Fig. [T] We see that, as photons diffuse over smaller 
distances, as a result of higher local baryon density in the 
presence of a CIP with positive A, the transition to ex- 
ponential damping of CMB anisotropies occurs at higher 
I. 

In the bottom panel of Fig. [U we show the visibil- 
ity functions g{z) = e~'^dT/dz for 3 different values of 
A; r is the optical depth due to Thomson scattering. 
The peak of the visibility function is the redshift zsls, at 
which most CMB photons last scatter. In the presence 
of a positive (negative) A CIP, decoupling occurs later 
(earlier) due to higher (lower) baryon density. 

To calculate the effects on the CMB moments Tim, we 
perturb the line-of-sight solutions, Eqs. dH and (|3]). This 
approach is relatively simple and amenable to rapid com- 
putation. The results should be accurate for multipole 
moments L < 870 for the CIP, as the baryon fluctuation 
can be considered as roughly constant in a given direction 
n across the thickness of the surface of last scatter on such 
scales. We discuss in Sec. IIV Al below how we extrapolate 
these results to smaller angular scales {L > 870) with a 
Limber approximation. 

We proceed by Taylor expanding in real space: 

f{v,n) = /(°)(77) + A(n)^(77) 



dA 



(8) 



where /^^'{t]) is the value of / under the null hypothesis 
A(n) — 0. We expand 



A(")"=H^i*^^^*^(")' 



(9) 



LM 



in terms of spherical-harmonic coefficients Alm for the 
angular variation in the CIP at the surface of last scatter. 
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FIG. 1: Dependence of physically relevant scales and Thomson scattering visibility function g[z) on amplitude of a global CIP 
perturbation A. Top left panel shows dependence of the angular sound horizon Is on A. Top right panel shows dependence of 
the diffusion damping scale Id on A. Bottom panel shows g[z) evaluated for 3 different values of A. 



We then apply the expansion in Eq. ([8]) to linear order in 
A to the line-of-sight expression, Eq. ([2]), and integrate 
over angles to obtain the first-order correction, 
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to Ti„i in the presence of a CIP, where 
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(11) 

(12) 



and the arrays inside parentheses are Wigner-3J sym- 
bols. Throughout, we use the indices L and M exclu- 



sively for the decomposition of the CIP, while lower-case 
indices are used for the multipole moments of the CMB 
observables Tim, Eim, and Bim- Sums over m (M) are 
always taken over the range —l<m<l {—L < M < L), 
while sums over / (L) are formally taken over 1 < Z < oo 
(1 < L < oo); in practice, a maximum value ^max is 
used for numerical evaluation, as discussed in Sec. IIVI 
The monopole L = corresponds to a global shift in fli, 
and ric and we absorb this term into the cosmological 
parameters themselves. 

For a given realization of the CIP field — that is, for a 
given set of A^jvf — the covariance between temperature 
moments is now 



(^Z'm'^im) — Ci 
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(13) 



of CIPs on the CMB power spectrum. We thus obtain, 
to second order in A, the temperature power spectrum. 
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where r/„j is the unperturbed moment in Eq. ([3]) , 5Ti^ is 
given by Eq. PH)) . and C; is the unperturbed power 

spectrum given by Eq. ([S]). The ^^^ superscript denotes 
the term arising when expanding Tim to order A^(fi). 
We evaluate this term using the second-derivative term 
d'^f^°'> /dA^ in Eq. (|C2p . We take an expectation value 
over CIPs and primordial-potential realizations. We then 
use Eqs. ©, dfl]), and ^7^, identities of Wigner-3J sym- 
bols [6g|, and Appendix [Cj to obtain 



Sy^^{cJ^'^ + Cr'^)Kt„ (14) 



■';;' 



and 



^X,dX' ^ 2 /■fc2^^p^^^-)^^(-^-) 



dA ' 



(15) 



for {X, X'} e {T, E, B} (a generalization that will be use- 
ful below), and C; is the temperature power spec- 
trum in the absence of CIPs. Here, 



— — = / dii—;^3i[Km - m 



dA 



dA 



(16) 



describes the change in the transfer function, Eq. ([6|, 
with A. 

In deriving these results, we have taken an average over 
realizations of the primordial-potential power spectrum 
$j:, but we have restricted our consideration to a given 
realization of the CIP. In Sec. |Vl we build the formalism 
to reconstruct A{n) from these off-diagonal temperature 
correlations as well their generalization to polarization. 
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The CIP power spectrum C^ and total variance A^^ 
obey 
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while the CMB derivative power spectra are given by 
.dx,dx' 2 /•,2„, „ „.dXi{k)dX[{k) 
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D. Temperature anisotropies with CIPs: Average 
over CIP realizations 



where d^Xi/d^A are defined analogously to the first- 
derivative transfer function in Eq. (|16p. Appendix |B] de- 
tails the calculation of the derivative power spectra. 



We now take an ensemble average over many realiza- 
tions of both the primordial potential field and the CIP 
field. This allows us, given a spectrum of CIPs, to calcu- 
late the effects of these CIPs on the power spectrum of 
CMB fluctuations measured on the entire sky. 

We denote the ensemble average of a spatially varying 
field X over realizations of the CIP field by {X)-^. We 
denote the ensemble average over realizations of both the 
CIP field and the primordial potential by (A) be- From 
Eqs. dnD and (m, we see that {T*^,TiJ\,^ oc (ALA/)b. 
For an isotropic random field, (AiM)b = 0, so we must 
thus go to second order in A{n) to compute the effects 



III. PERTURBED LINE-OF-SIGHT 
FORMALISM: POLARIZATION 

We now generalize the analysis above to the CMB po- 
larization. In addition to inducing off-diagonal correla- 
tions in the polarization spherical-harmonic coefficients, 
CIPs will induce B modes. We begin by reviewing the 
LOS solution for polarization under the null hypothesis of 
no CIPs. We then compute the effects of CIPs, both for 
a single realization of the CIPs, and then for an average 
over realizations of a spectrum of CIPs. 



A. Polarization anisotropies with homogeneous 
baryon fraction 

Polarization is a spin-2 tensor field and can be ex- 
panded as [63, [el] 



Pab (n) = Y, {EimYil^,, + Bi^Y^^^J , 



(24) 



Just as in the case of temperature, when considering a 
single realization, we need only consider the first-order 
terms in Eq. ([29]). We then utilize Eqs. ^, ^, and 
the first-derivative piece of the usual expansion for A.{h) 
[see Eq. ([8])] to obtain an expansion for the polarization 
tensor Pai in the presence of CIPs: 



Irn 



where YJ^ ^^ and Y;^ ^^ are the E- and B mode ( "grad" 
and "curl" , respectively) tensor spherical harmonics, as 
defined in Appendix El The right-most indices after the 
comma are tensor indices. In terms of the Stokes polar- 
ization parameters Q and U, the polarization tensor is 
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Pab (n) = I 



Q [n] — t/(n)sin6 

-U{n)s\n9 — Q (ri) sin^ ( 



(25) 



where 9 is the polar angle of the LOS with respect to 
some origin. 

Under the null hypothesis, the polarization pattern at 
the surface of last scatter is a pure E mode, with multi- 
pole moments given by 



We may now pick off the induced E- and B-modc mul- 
tipole moments SElJ^ and 5B\2_ at order A, using the 
appropriate integral over a tensor spherical harmonic: 



47ri' /■ f 

Elm = —-3 j d^ f^{^) J d^k^j:ji{x)Yi*„Ak), (26) 

where x = fc(?yo ~ ^7), and f^{r}) is the E-mode transfer 
function, obtainable numerically from Boltzmann codes. 
The polarization covariance and TE covariance are de- 
rived analogously to the results for temperature, yielding 
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We evaluate Eqs. ([32]) -((Ml), calhng on Eqs. (IST0l)-(lM2|, 
yielding 
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B. Polarization anisotropies w^ith CIPs: single CIP 
realization 



where 



We now generalize the analysis to include the effects of 
CIPs. In the presence of a CIP field A{ii), the real-space 
polarization tensor may be Taylor expanded as 
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We now evaluate the induced correlations between dif- 
ferent temperature/polarization moments. At first order 



■mA{n),{Bp^,Bi^)^{B;[ll6Bl'J) 



0. The remaining 



covariances are 



power spectra are 



L-\-l-\-l' even 



LM cL,EE 



LM 
r,L,EE _ /^E.dE , ^E,dE\ ttL 



'W ' 



L+l+l' odd 

/P* R \ — Y^ A C^^^t qL,EB 

K^l'm'^lm) — 2^ ^LKKlm^l'm'^lV ' 

LM 



L+i+i' odd 



(37) 



(38) 



LM 

Slf'^^-^C^''^Ht,. (39) 



L+l+l' even 



EA f:LM qL, 

^LMC.lmMm'^W 



X,TE 



LM 

■~1 



^.^;™- {cr''Hli,+Cf'^-Kt,). (40) 



C. Polarization anisotropies with CIPs: Average 
over CIP realizations 



We now extend the ensemble average to multiple real- 
izations of the CIP field. We do this in order to compare 
the polarization power spectrum induced by CIPs at the 
surface of last scatter with that induced at reionization. 
For temperature, the average over realizations of both the 
CIP and primordial-potential perturbations is given by 
Eq. (dZl). Extending this average to X,X' e {T,E,B}, 
we obtain the XX' power spectra, averaged over the en- 
tire sky, to second order in A: 
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\ /be 
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sxir^i 



L+li+l even 



■"l 



^l 



(42) 



^^TE,(l) 



'■ 2 



E,-rA,-rdT,dE^ rjL jzL 
L,h 



[cj-"'"' + c^-"'^ 



Cr ^ C™'("' + SCf^-^'^ + JC™'('\ (43) 

L-{-li-\-l even 

^^EE,i)^ ^ Ctct^-'''GuAHii)\ 

LJi 



JC™'(2' = Alcf^"'-^, 



and 



L+h+l odd 



Cf^c. Y. Ctcf^'^^GuAHliy ■ (44) 

Ldi 

The CIP field A(n) is a scalar and cannot statistically 
change the parity of polarization perturbations. This re- 
quires that C^^ and C™ vanish when averaging over 
CIP realizations. Algebraically, this is enforced by the 
vanishing of the relevant Wigner-3 J symbols, as occurs 
with optical-depth fluctuations at reionization [44l . |69| 
and wit h g ravitational-potential perturbations in weak 
lensing (toTItiI 



Indeed, the geometric (Wigncr-3 J) sym- 
bols obtained are the same as for those effects. CIPs give 
rise to different W dependences for the functions S*;;,' , 
however, through the dependence on the derivative power 
spectra Cf^-^ , allowing them to be disentangled obser- 
vationally from gravitational-potential fluctuations along 
the LOS or optical-depth fluctuations at reionization. 



IV. 



NUMERICAL RESULTS FOR B- MODE 
POWER SPECTRA 



We now apply the formulas derived in Sees. IIIDI and 
unci to compute the power spectra for B modes induced 
by CIPs at decoupling. We first discuss the form of 
the angular CIP power spectrum C^. We then present 
numerical results for B modes induced at decoupling. 
For comparison, we then reproduce the calculations of 
Ref. 133 of the B modes induced at reionization. 



where C, 



XX', (0) 



is the power spectrum computed with no 



CIP contribution. 



We evaluate Eq. (gl]) with Eqs. (|A12|) and (|A13|) and 
Wigner-3J relations to simplify the resulting integrals 
and sums. Superscript indices ^^' and ^'^' indicate the 
order of the derivative cPf'-^^'/dA^ used to derive the 
indicated term, as in Sec. IIIDI The resulting nonzero 



A. Power spectrum of compensated perturbations 

1. Three-dimensional CIP power spectrum 

To proceed further, wc must make an ansatz for the 
spectrum of CIPs. Motivated by the curvaton model 
(which produces a nearly scale-invariant spectrum of 



CIPs) [SJ, [Sa, E^l 1 we assume a seale- invariant spectrum 
for the three-dimensional CIP field A(x); that is, 



(45) 



A*(k')A(k)^ M^nfSl (k - k') PA(fc) 
PA{k) - Ak-\ 



where A(k) is the Fourier transform of A(x), and A is a 
dimensionless CIP amplitude. 

As discussed in the introduction, the strongest con- 
straint to A comes from cluster baryon fractions. This 
constraint tells us that the variance, 






k'^dk[3ji{kR)/{kR)fPA{k), (46) 



in the baryon-to-dark-matter ratio on i? ~ 10 Mpc scales 
is Aci < 0.08. The integral has a formal logarithmic di- 
vergence at low k which is cut off, however, by the vol- 
ume occupied by the clusters surveyed. Taking this to 
be the horizon, fcmin — (10 Gpc)^^, we find A < 0.017. 
Since the cosmological baryon fraction fl^, determines pri- 
mordial abundances via BBN, there is an additional con- 
straint from astrophysical measurements of these abun- 
dances [33 • However, this constraint is less stringent 
than the one from cluster gas fractions. 



computation in Ref. [7^). Moreover, as we shall see 
in Sec. I VII most signal-to-noise in CIP reconstruction 
comes either from L < 100 or L ^ 2000, and so the main 
conclusions of this work should not be affected. 
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2. Angular CIP power spectrum 



When the 3-dimensional field is projected onto a nar- 
row spherical surface, the resulting angular power spec- 
trum for A will be C^ ~ A/{t:L'^) for mulipole moments 
-^ ^ (% ~ ^is)/o'r; — 870, where r/ig and 770 are the con- 
formal time at last scatter and today, respectively, and 
Gn is the rms conformal-time width of the surface of last 
scattering (SLS). At smaller angular scales (larger i), 
the angular variation in A is suppressed by the finite 
width [74I of the scattering surface. Using the Limber 
approximation, the angular power spectrum for A can 
be approximated by C^ ~ ^(??o — Vis)/{'2\/TTL^cr,j) for 
L > 870. The exact analytic expression we use is ob- 
tained from the Limber approximation, approximating 
the visibility function as a Gaussian. It is 






u 



2 \?7o 



La„ 



r\s\s 



(47) 



where t/(a, 6, x) is a confluent hypergeometric function. 
We use 770 — 771s = 14100 Mpc and a^^ — 16.2 Mpc for de- 
coupling. We use ?7o — r\\s — 9760 Mpc and a^^ = 448 Mpc 
for reionization. These values are obtained by directly fit- 
ting to the visibility function output by the CAMB code 



6l|. Of course, Eq. (|T7)) is an approximation, and the 



precise shape of the transition from C^ ex Xjl? — > Cf 
1/L^ near L ~ 870 depends on the interference of Fourier 
modes of A with those of $, averaged over the SLS. This 
issue warrants future study, but the asymptotic behavior 
at low and high L is correct (as shown for an analogous 



FIG. 2: CMB B-mode polarization power spectra induced 
by CIP perturbations at decoupling (black solid line), com- 
pared with the effects of CIPs at reionization, for which two 
contributions are shown: patchy screening (red dotted line), 
and patchy scattering (blue short-dashed line). The ampli- 
tude for the CIP power spectrum is that which saturates the 
Aci < 0.08 bound from clusters [3q |. 



B. Numerical result for B modes from CIPs at 
decoupling 

Using the Limber approximation with values 770 ~ '7is = 
14100 Mpc and cr,, = 16.2 Mpc appropriate for decou- 
phng, Eqs. (|42|) ~ (|44| can be used to obtain predictions 
for the B modes induced by CIPs at decoupling. The 
results arc shown in Fig. [2] for A = 0.017, the largest 
CIP amplitude consistent with the galaxy-cluster con- 
straint. Appendix IBJ details the calculation of the requi- 
site derivative power spectra. We use a maxinrum I value 
of Z„iax = 10000. 



C. Reionization 

In Ref. [33, it was noted that spatial inhomogeneities 
in the baryon density give rise to angular variations in 
the optical depth r for rescattering of CMB photons 
at reionization. It was also noted that these inhoino- 
geneities would give rise to B modes primarily at large 
angular scales by patchy rescattering of CMB photons 



and at smaller angular scales through patchy screen- 
ing of the primary CMB polarization from the decou- 
pling epoch. These calculations build upon calculations 
in Refs. [4J, |6^, |74h76| where optical-depth fluctuations 
were postulated to arise from inhomogeneities in the free- 
electron fraction due to inhomogeneous reionization. 

In our notation, the contribution of patchy screening 
is H i69,l 






(48) 



L-\-li-\-l even 



(49) 



L-\-l\+l even 
jrriEE 2 -2t Y^ ^A/-.EE,rcc /tjL \2 „ 

L.h 

(50) 

L+h+l odd 

SCr =r^e-^^ Yl Ctcf^^- «)' Gu., (51) 

where r = 0.086 is the mean optical depth, and Cf' 
is here the angular CIP power spectrum for reioniza- 
tion; i.e., obtained with 770 — rjis — 9760 Mpc and 
ari = 448 Mpc. These values are obtained by fitting a 
Gaussian visibility function to the reionization model of 
Ref. [73. 

The contributions of patchy scattering are 



SC, 



BB 



e^EE _ •J" ,^A,o2 

^^' ^100^' 



•'/ 



--/ 



(52) 



where Qrms — 17.9 /iK is the rms temperature 
quadrupole at reionization. 

Figure [5] shows the B modes induced by patchy scatter- 
ing and screening at reionization again using A = 0.017. 
We see that at all but the largest scales, the decoupling- 
induced B modes are larger (by up to ^ 3 orders of mag- 
nitude) than those induced at reionization. 

We thus conclude that the effects of CIPs on CMB fluc- 
tuations would be much larger than found in Ref. [35|, 
particularly at the small scales most important for detec- 
tion and reconstruction (to be discussed below) of CIPs 
from the CMB. We thus now move on to show how spa- 
tial fluctuations in the baryon-to-dark-matter ratio can 
be measured with CMB maps. 



V. MEASUREMENT OF CIPS WITH THE CMB 

In this section, we show how the CIP field A{n) can be 
measured with off-diagonal CMB correlations, building 
upon analogous prior work on measurement of cosmic- 
shear fields iTOl. [71I [78l - [8l| . departures from statisti- 
cal isotropy [8a - l84| . and cosmic birefringence [85l - l88| . 
Having concluded that the decoupling signal is much 



TABLE I: The "response functions" Sf;;^^' of CMB fluctu- 
ations to CIPs, defined in Eqs. ITS]) and ([371) -(gO]), for the 
various correlation functions. 



XX' 



qL,XX' 



TT 
EE 

EB 

TB 

TE 



T.dT^^T.dT^)^.^^^ 



^E,dE . ^E,dE\ rjL 
C; + C,, j U^, 

-jG,, Nil, 



C, 



T,dE^.^^E,dT^. 



bigger than that from reionization, we consider detec- 
tion/measurement of CIPs at the surface of last scatter. 

In Sec. IV Al we construct a minimum- variance 
quadratic estimator A^m for the multipole moments of 
the CIP field. In Sec. IV Bl we explicitly calculate the 
noise covariance (due both to cosmic variance and exper- 
imental noise) needed to evaluate the errors and optimal 
weights of Sec. IV Al Finally, in Sec. IV CI we use the re- 
sults of the preceding sections to derive an expression 
for the signal-to-noise ratio (SNR) with which a CMB 
experiment can detect CIPs. 



A. Minimum-variance estimators for Alm 

The total correlation between multipole moments (in- 
cluding the contribution induced by a given realization 
of CIPs) takes the form [see Eq. ^ and (P7 1) -(|ig i) ]. 
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Ci Sw 6mm' 
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(53) 



LM 



where 



^LM,XX' 

hi' 



^L.XX' 

hi' 



As before, X,X' G {T,E,B}. As discussed in Sees. Ill 
and mil the functions S^! map A^a/ to observed off- 

diagonal CMB anisotropies. The S^, ' are assembled 
in Table H 

The spherical-harmonic coefficients X"^'^^ obtained by 
a given CMB experiment are related to the true coeffi- 
cients Xim by Xj'^fP = WiXim, where Wi = e-'('+i)-'/2 
is the window function, and ai, ~ ^fwhm/V8 In 2 ~ 
0.00741 (6'f„hm/l°) is related to the beam's full width at 
half maximum (FWHM) 6'fwhm- The observed two-point 
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correlations are then 
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(55) 



Following Refs. |44l. I69l. I84| - l86l l88j . the minimum- variance 

quadratic estimator for the rotational invariant D^, ^^ 
is 






LA/ 
Ira.V 



(56) 



is odd. The appropriately normalized covariance matrix 

for 5i^^^AA\map j^ 



CAA' /-» / / p»LA/.A,niap ^LA/,A'.rQap* 



(59) 



^LJ\f, A. inap\ / p;LA/.A',map* 



We now proceed to compute the covariance matrix 

/^AA' 



B. OfF-diagonal covariances 



To extract Ala/, CMB temperature and polarization 
maps must be used to reconstruct 0^, ' '™^^ by ap- 
plying Eq. ([56|) . Then, through the estimator A^j{,^ = 
j5LM,xx',map^/^^^^^^L,xx'\^ we obtain many mea- 
surements of Aij\/. These measurements are generally 
correlated (even for fixed UV). so we must take care to 
construct an optimal estimator A^a/ for CIPs when using 
the full set of available maps for a given experiment. Gen- 
eralizing the estimator and error formulae in Refs. J84ll8a | 
to our case of interest, we obtain the optimal estimator 
AiM and its rms error a/^^ , taking into account all pos- 
sible correlations between X and X': 



To numerically evaluate Eq. (|59p . we must have a 
model for the statistics of the observed map covariances, 
including noise. We assume the noise in each pixel is sta- 
tistically independent, Gaussian, and uncorrelated with 
the signal, and we assume no coupling between the noises 
in {T, E, B}. In this case, the noise power spectra are 



^BB, noise ^EE, noise 



= 2C, 



TT. noise 
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Ala/ = cxi, Y. ^iv E ^u'^ 5^ ■^^■"^ {C^A 

l'>l AA' 

a^l = E Gil' E ^li'^ ^li'^ ['^u'^] AA' ' (5^) fo'' tli° ™aP are 



where NET is the (effective) noise-equivalent tempera- 
ture for the experiment, iobs the duration of the exper- 
iment, and /survey IS the fraction of sky surveyed. The 
cross spectra C; ^"°'^° = if X 7^ X' . The power spectra 



i'>i 



AA' 



where {A, A'} e {EB,BE,TB,BT,TT,EE,TE,ET} 
when I ^ I, {A, A'} e {EB, TB, TT, EE, TE} when I = I', 
and 
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The inequality I' > I is imposed so that we do not double 
count correlations. Sums are subject to the additional 
restriction that for {A,A'} € (TE,ET, EE,TT), l + l' + L 
is even, while for {A, A'} e (BE, EB, BT, TB), l + l' + L 



is useful to explicitly evaluate Eq. (|59|) using 
Eq. (|56|) . Since all fields involved are Gaussian, all the 
four-point functions that arise may be evaluated using 
Wick's theorem. Wigner-3 J identities may then be fruit- 
fully applied to obtain all the elements of Cj'^A expressed 
in terms of C; '™3,p^ j£ ^ _ j^/ ^ ^i^^j^ C;z is a 5 x 5 diago- 
nal matrix, with rows/columns in the order TT, EE, TE, 
BE, BT, and entries 
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If / 7^ I' , Cw is an 8 X 8 block-diagonal matrix, with 
rows/columns in the order TT, EE, TE, ET, BE, EB, 



BT, TB, and entries 



11 



Mw 







K-w 



Ml 



/ ^TT,map^TT,map ^TE,map^TE,map ^TT,map^TE,map ^TE,map^TT,map \ 

^TE,map^TE.map ^EE,map^EE.niap ^TE.map^EE,map ^EE,map^TE.map 

^l '-'l' '-■l '-'l' '-^i '-'I' '-■l '^l' 

^TT,map^TE,map ^TE,map^EE,map ^TT,map^EE,map ^TE,map^TE,inap 

•-"; '-^v ^i ^i' <-'; <-';' ^i '-^i' 

^TE,map^TT,map ^EEjinap^TE.niap ^TE.inap^TE,map ^EE,map^TT.inap 

\ '-'I '-'V ^i '-'V <-^i '-^v '-^l '-'V 



, (64) 



/ ^BB.map^EE 



/C/c = 



^BB , map ^TE , map 



3^T 










^EE 


map^BB 


map 









c™ 


map y-yBB,m 


^BB 


map^TE 


map 









cr 


map^TT 


m.ap 














^TE 


inap^BB 


map 









y^TT/^BB.map 



r 



In Sec. lVI Al we apply the preceding formulae to estimate 
the noise in the reconstructed CIP field for a variety of 
ongoing and upcoming experiments. 

In many cases, most of the sensivity to CIPs comes 
from a single combination (e.g., TT or TB) of observ- 
ables. It is therefore interesting to consider the constrain- 
ing power of a single such combination. In the case of TB, 
the error is given by [see Eqs. ()57p ] 



l+l'+L odd 

E 
i'>i 



Gw 



S, 
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(65) 



To generate the noise curves discussed in Sec. IVI A| we 
use expressions analogous to Eq. (|55|) for each pair of ob- 
servables. These noise values are then applied to estimate 
the SNR with which a given spectrum of CIPs might be 
detected. 



C. Signal-to-noise formula 

Ultimately, we wish to assess the SNR of our estima- 
tors for a given CIP power-spectrum amplitude A. Each 
estimator A^m gives an independent estimator for A, 
and by adding them all with inverse-variance weighting, 
the total SNR with which CIPs can be detected is 
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/sky 
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(66) 



where the error is evaluated using Eq. ([57)1 . /sky is the 
fraction of sky used in the data analysis, and Lmin = 

— 1/2 

/sky • Modes that vary on scales larger than the area 
of sky analyzed will have degraded signal-to-noise. A 
minimum value of L is thus imposed to conservatively 
account for fractional sky coverage. In practice, these 
modes would still contribute to the integrated CIP power 
in the area of sky analyzed. In this work, however, we im- 
pose a cut at Lmin to avoid an over-estimate of sensitivity, 
all the same establishing the utility of CMB observations 
for probing CIPs. 



VI. EXPERIMENTAL PROSPECTS 

We now apply the formalism of Sec. |V] to assess the 
prospects of using CMB experiments to detect CIPs. We 
consider specifically the ongoing WMAP [90] and Planck 
[43 1 satellites and a possible future satellite, EPIC [46| . 
We also consider the following suborbital experiments: 
Polarbear ^, Spider ^ and ACT [gland SPT ^ 
and their polarization upgrades ACTPol [4j] and SPTPol 
[45| . Finally we consider an idealistic cosmic- variance 
limited (CVL) experiment, limited only by sky cuts to 
avoid galactic foreground emission. We do this to roughly 
quantify the lowest CIP amplitudes that could ever be 
probed with the CMB. 

The experimental parameters assumed for these exper- 
iments are given in Table II. For WMAP, Planck, EPIC, 

= 1, while 



and the CVL case, we assume that /^ 



survey 



for Polarbear, SPT, ACT, SPTPol, and ACTPol, we as- 
sume that /survey = /sky For WMAP, we assume use 



12 



of the V and W bands in the analysis. For Planck, we 
assume that the 143 and 217 Ghz channels will be usable 
for CMB anisotropy measurements, and take appropriate 
inverse- variance weighted sums of the noise in these chan- 
nels. Appendix |B] details the calculation of the requisite 
derivative power spectra. We include BB correlations 
induced by gravitational Icnsing when evaluating C^,^ , 
using the CAMB Icnsing module J6l| . 
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WMAP 


W Band 
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0.65 
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Polarbear 
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Planck HFI 
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0.65 


1.2 


Planck HFI 


217 Ghz 
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91 


0.65 


1.2 


Spider 


150 Ghz 


30 


4.2 
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0.016 


ACT 


148 Ghz 


1.4 


58 


0.0072 


0.14 


SPT 


150 Ghz 


1.2 


91 


0.0024 


0.29 


ACTPol 


150 Ghz 


1.4 


6.0 


0.10 


0.21 


SPTPol 


150 Ghz 


1.0 


14 


0.016 


0.75 


EPIC 


150 Ghz 


5.0 


2.0 


0.65 


4.0 


CVL 






0.0 


0.65 





TABLE II: Experimental parameters for the experiments con- 
sidered in this work: beamwidth 9 (in arcminutes), noise- 
equivalent temperature (NET) (in fiK sec^"), and observa- 
tion time fobs (in years). 



A. Noise curves 

We compute the noise curves in the reconstruction 
^LM for all experiments under consideration, and show 
the results in Figs. [SHSl Wc plot the noise power spec- 
trum 



SCl 



fs^yV2LTT 



(67) 



as well as the signal Cl- We use the value A = 0.017, 
which saturates the galaxy-cluster bound on CIPs. Ex- 
periments with larger beams, such as Polarbear and Spi- 
der, generally have higher noise levels (for reconstruction 
of A) than do others. At low L, temperature is generally 
the best probe of CIPs. At sufhciently high L, the BT 
correlation takes over as the best probe of CIPs. At very 
high L, all the noise curves grow very large, indicating 
that the fidelity of the reconstruction of the CIP breaks 
down at small scales. This is very similar to lensing and 
is expected, as fluctuations in the baryon fraction couple 
small to large scales. We note that C^ ''"^^ is computed 
using T — 0.086, as the observed CMB anisotropics are 
affected both by a screening envelope at high I, due to 



reionization, and by a reionization "bump" in polariza- 
tion at low /. We use CAMb's built-in tank reionization 
model with the parameters of Ref. [73] . 

We then compute the total noise in the reconstruc- 
tion A^M for each experiment, adding the different cor- 
relations in quadrature with inverse-variance weighting. 
This should be a reasonable approximation to the sum in 
Eq. (j57p , as inverse- variance weighting tends to be domi- 
nated sharply by the lowest noise correlation. The results 
are shown in the left panel of Fig. [71 We sec that broadly 
speaking (with some alternation as a function of L), the 
best sensitivity is achieved by EPIC, followed by SPT- 
Pol, ACTPol, Planck, Polarbear, WMAP, ACT, Spider, 
and SPT. Obviously, any specific experimental concept 
will be outperformed by the CVL case, as confirmed in 
the left panel of Fig. [71 



B. Signal-to- noise 

Calling on Eq. (|55)) . we compute the SNR expected for 
all the CMB experiments wc consider as a function of the 
rms CIP fluctuation A^. The results are shown in the 
right panel of Fig. [71 Assuming a scale-invariant spec- 
trum, we see that already with WMAP, we are able to 
probe compensated fluctuations in the b ary on fraction of 
^ 10%, the lowest value probed in ReL [Sa]- We dcflne a 
"detection" as a measurement with S/N > 3. Currently 
operating suborbital experiments, like SPT and ACT, 
and the upcoming Spider experiment, perform compara- 
bly to WMAP. Although these experiments are sensitive 
to a fairly large rms CIP amplitude, the cluster constraint 
was obtained at a different scale, and it is important to 
check the constraint using the independent probe offered 
by the CMB. 

Planck and Polarbear offer the next major improve- 
ment in SNR, probing Ad ~ 3% and higher. The addi- 
tion of polarization sensitivity to SPT and ACT lowers 
the range of detectable Ad to 0.6%. EPIC would be able 
to measure Ad ~ 0.4%. This is a factor of 20 lower than 
the currently allowed maximum from measurements of 
the baryon fraction in galaxy clusters. In the CVL case, 
an additional order of magnitude improvement in SNR is 
possible. 

As discussed in Sec. IV Cl we conservatively estimated 
sensitivity, omitting low-L modes. In practice, all-sky 
experiments like WMAP, Planck, and EPIC could probe 
even smaller CIP amplitudes than estimated here, as the 
L = 1 mode could contribute signiflcantly to the SNR. 
If experimental techniques improve (approaching the 
cosmic- variance limit at very high I) and a way is found 
to disentangle CIPs from secondary CMB anisotropics at 
these high /, a further order of magnitude improvement 
in sensitivity is theoretically possible, using the estima- 
tors developed in this work. Additionally, we have been 
conservative in our estimates of SNR, assuming that only 
one useful frequency channel is available for EPIC. The 
mission concept actually calls for ^ 7 channels, in order 
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FIG. 3: Predicted noise power spectrum &Cl in the reconstructed CIP perturbation field A_la/ in four different ongoing/future 
CMB anisotropy experiments, as a function of angular scale L. We separately plot the noise for distinct pairs of observables: 
TT is shown as an orange (short-dashed long-dashed), TE eis a green (dotted-dashed) line, EE as a magenta (short-dashed) 
line, BE as a blue (long-dashed) line, and BT as a red (dotted) line. Also shown (black solid line) is the power spectrum C^ 
, marked signal, for a scale-invariant spectrum of CIPs with the maximum amplitude allowed by galaxy clusters. Each panel 
shows estimates for a different experiment, as indicated in the figure. The beige (grey) shaded region shows the range of L that 
is not included in our estimates, due to finite sky coverage effects, as discussed in Sec. IV CI 



to achieve good foreground subtraction. It may be that 
if EPIC is built, more than one useful channel of signal is 
obtained, improving the SNR by a factor of order unity. 
It may be that CIPs are not an independent Gaussian 
random field (as assumed here), but rather, as in some 
curvaton models fSJl , correlated with the usual adiabatic 



fluctuations. In that case, CIPs will induce 3-point cor- 
relations between CMB observables, through effects at 
decoupling. In future work, we will pursue the possibility 
of probing CIPs with the corresponding CMB bispectra. 
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FIG. 4: Predicted noise power spectrum SCl in the reconstructed CIP perturbation field Alm in four difi^erent ongoing/future 
CMB anisotropy experiments, as a function of angular scale L. Colors (line styles) are as in Fig. [3] We separately plot the 
noise for distinct pairs of observables: TT is shown as an orange (short-dashed long-dashed), TE as a green (dotted-dashed) 
line, EE as a magenta (short-dashed) line, BE as a blue (long-dashed) line, and BT as a red (dotted) line. Also shown (black 
solid line) is the power spectrum C^ for a scale-invariant spectrum of CIPs with the maximum amplitude allowed by galaxy 
clusters. Each panel shows estimates for a different experiment, as indicated in the figure. The beige (grey) shaded region 
shows the range of L that is not included in our estimates, due to finite sky coverage effects, as discussed in Sec. IV CI 



VII. CONCLUSIONS 



Compensated isocurvature perturbations provide an 
intriguing empirical possibility for large-amplitude de- 
partures from homogeneity in the early Universe. The 
current constraint, < 10%, to the amplitude of such 
perturbations is surprisingly weak. While Ref. |35| has 
pointed out that there may be CMB signatures induced 
by CIPs at reionization, we have shown here that the 



CMB effects of CIPs at the surface of last scatter would 
be several orders of magnitude larger. We then calcu- 
lated the full two-point temperature/polarization corre- 
lations induced by CIPs on the CMB and developed the 
minimum- variance estimators for measuring the CIP field 
with the CMB. 

The WMAP satellite may be sensitive to a scale- 
invariant spectrum of CIPs, but only if the CIP ampli- 
tude is close to its current upper bound. In the future, 
sensitivity to CIP amplitudes as small as ^ 3% may be 
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FIG. 6: Predicted noise power spectrum SCl in the recon- 
structed CIP perturbation field Alm for an ideal cosmic vari- 
ance limited experiment, as a function of angular scale L. 
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lowed by galaxy clusters. The beige (grey) shaded region 
shows the range of L that is not included in our estimates, 
due to finite sky coverage effects, as discussed in Sec. IV CI 



achieved by instruments in operation and ^ 0.1%-level 
fluctuations accessible in the near future with precise 
ground- and space-based polarization experiments that 
are under construction (ACTPol and SPTPol) or in con- 
ceptual development (EPIC). 

Many steps must be taken before such measurements 
can be implemented with real data. Techniques to 
deal with partial-sky coverage and realistic instrumental 
noise properties must be developed, but these techniques 
should be similar to those being developed already to 
measure the effects of weak gravitational lensing on the 
CMB. Likewise, techniques must be developed to distin- 
guish the off-diagonal correlations induced by CIPs from 
those induced by weak gravitational lensing (e.g., j95l|). 
which should be comparable in amplitude if the CIP field 
A is comparable to the lensing potential (j), i.e, ^ 1%. 

Although lensing is already included in our reconstruc- 
tion noise estimates [ctal in Eq. ([ST]) ], it might also 
induce bias in the reconstruction of the CIPs. Gravi- 
tational lensing of the CMB will induce correlations of 
similar form to Eqs. ([55)1 and (|54l) . with the distinction 

that the functions 5;;/ will be different for lensing than 
for CIPs. In the case of lensing, these functions describe 
the remapping of CMB observables on a lensing-deflected 



sky. In the case of CIPs, these functions describe the de- 
tailed physical dependence of CMB anisotropics on the 
baryon density. If Eq. ([57| is applied to the extension of 
Eqs. ([53)) -([54 l) that includes gravitational lensing, a bias 
will be induced in the measurement of Alm- 

The differing forms of S*;;,' will allow lensing and 
CIPs to be disentangled. Using a straightforward gener- 
alization of the estimator in Eq. ([57)) that includes terms 
for both CIPs and lensing, the bias induced by lensing 
on CIP measurements may be removed, simultaneously 
reconstructing the lensing and CIP fields. This is anal- 
ogous to the estimators discussed in Ref. [95[, where it 
is shown that if reionization is patchy because of inho- 
mogeneity in the distribution of ionized bubbles around 
the first sources, the contributions of patchy reionization 
and lensing may be distinguished. In future work, we will 
explicitly compute the bias in CIP searches that will be 
induced by lensing and will construct the estimators that 



disentangle lensing from CIPs. As in Ref. [95[, we expect 
that the signal-to-noise of the biased and unbiased esti- 
mators should be nearly the same, and so gravitational 
lensing should not affect the signal-to-noise estimates of 
this paper. 

The measurements we propose in this paper offer a 
precise test of how closely the primordial baryon and 
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dark-matter distributions are matched and approach the 
CIP amplitudes allowed in curvaton models. Moreover, 
if future CMB experiments detect subdominant isocur- 
vature fluctuations between matter and radiation, the 
techniques developed in this work could disentangle con- 
tributions from the baryon and cold dark-matter (CDM) 
isocurvature modes. Even greater gains in sensitivity are 
theoretically possible if the cosmic- variance limit is ap- 
proached at high I by future experiments. We are op- 
timistic that, in the near future, we will learn just how 
well baryons trace dark matter in the early Universe. 
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Appendix A: Tensors on the sphere 

Here we review the tensor spherical-harmonic formal- 
ism, following closely Ref. [8g|. The metric on the 2- 
sphere is 
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where d is the polar angle defined with respect to the 
origin of the orientation vector h. It is useful to introduce 
the orthonormal basis 







The tensor spherical harmonics are [67 
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where the normalization constant is given by 



m 



^ l 2{l-2)\ 

" V G + 2)! ' 



(A5) 



and all indices following ":" denote covariant derivatives 
taken on the 2-sphere. The indices I and m are multipole 
indices, while a and b are tensor indices. 

Using Ref. [701, the covariant derivatives may be ex- 
pressed in terms of spin- 2 spherical harmonics |68l. uOj : 
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The left subscript "2" denotes a spin-weighted spherical 
harmonic gY/m of spin s, while ® denotes a tensor prod- 
uct. The spherical basis vectors m± are 
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In row-column form, the spherical tensor basis functions 
are then [83 
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where wc have suppressed the Im indices for the sake of 
brevity. 

To evaluate the polarization anistropies induced by 
CIPs (see Sec. IIIip . it is useful to obtain identities for 
the product of two distinct (generally different I and m 
values) spherical harmonics. Using Eqs. (jA8p - (jA9p . it 
can be shown that (sdl 



(^E...)*yE 



(^B,a.)*yE 



{+2X* ® +2Y + _2^* ® ~2Y) , (AlO) 



= -- [+2X* ® +2Y - -2X* ® -2Y) , (All) 

where X and Y signify tensor spherical harmonics on 
the left-hand side of Eqs. (jAlOp and (jAlip and the corre- 
sponding spin- weighted spherical harmonics on the right- 
hand side of these equations. Relations such as these help 



express the integrals of Sec. IIIII as integrals over 3 spin- 
weighted spherical harmonics. These are then evaluated 
by applying the relation [70| 



(A12) 



In evaluating the contribution of terms proportional to 
(pf(^')/dA'^ to perturbed LOS solutions for polarization 
anisotropics, integrals over 4 spin- weighted spherical har- 
monics must be evaluated. They may be simplified using 
the identity ^ 
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Appendix B: Derivative power spectra 

Wc wish to estimate the derivatives cPXi (fc) /dA" of 
the transfer functions Xi (k). For the first derivatives, 
we use a 5-point approximation, running CAMB |61| re- 
peatedly to obtain the transfer functions at 5 different 
values of ^bh^ and rich^. The first derivative is then well 
approximated by 



dXi (fc) 
dA 



2 

E 

i=-2 



C^Xi (fc) 

12A 



(Bl) 



Here X^ (fc) denotes the transfer function evaluated under 
the transformation fib -^ 17;, (1 + A), fie -^ ^c ^ Afib. 
We find that the choice A = 0.02 works well to eval- 
uate the first derivatives. We run convergence tests by 
doubling and halving A and find that dXi (fc) /dA has 
converged to ^ 5%, which is more than sufficient for our 
purposes. We use values cq — 0, c±i — ±8, and c±2 = =Fl 

For the second derivatives, we use the 7-point numeri- 
cal approximation 



d^Xi (fc) 
dA2 



c^Xi (fc) 



'^ 180A2 



(B2) 



In this case, we find that the choice A — 0.066 lies com- 
fortably in the zone of convergence. The corresponding 
coefficients are cq = —490, c±i = 270, c±2 = —27, and 
c±3 = 2. We run convergence tests by doubling and 
halving A and find that d^Xi (fc) /dA^ has converged to 
^ 5%, which is accurate enough for the work presented 
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here. The resulting derivative power spectra, defined by 
Eqs. (I221)-(ES1), are shown in Figs. Q-TOl All derivative 
power spectra are computed using CAMB, with t ~ 0. 
These are then multiplied by a homogeneous reionization 
damping envelope with mean optical depth t ~ 0.086, 
given by expressions in Ref . |9M| • This was done to isolate 
the effects of patchy decoupling, screened by a homoge- 
neous optical depth at reionization (with ^reion — 10.5), 
from additional (smaller) anisotropics induced at reion- 
ization. 



Appendix C: Harmonic expansion of CMB transfer 
functions 



The most convenient way to generalize Eq. ^ to in- 
clude terms oc A^ (fi) is to derive second-order corrections 
to 



Using the Taylor expansion in real space defined by 
Eq. (HJ and Eq. (fCTj) . we obtain 
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